;***************************************************************
;discription：低温胁迫因子
;
;param1: nrows(input) 行
;param2: ncols(input) 列
;param3: t_opt(input)植被生长的最适宜温度
;param4: t_mon(input) 某一月的平均温度
;param5: e_t1(output) 低温胁迫因子
;
;author: wangqiang( wensten )  2015-09-14  chengdu
;***************************************************************
PRO PRO_E_T1, nrows, ncols, t_opt,t_mon, e_t1=e_t1
  
  ;开始计算每一个像元对应的低温胁迫因子
  FOR yindex=0,nrows-1 DO BEGIN
    FOR xindex=0,ncols-1 DO BEGIN         
      ;最适宜温度
      tmp_opt_cellvalue=t_opt[yindex,xindex]  
      ;月平均气温   
      tmp_mon_cellvalue=t_mon[yindex,xindex]
      
      IF (tmp_opt_cellvalue EQ -9999.0)||(tmp_mon_cellvalue EQ -9999.0) THEN BEGIN
        e_t1[yindex,xindex] = -9999.0
      ENDIF ELSE BEGIN        
        ;当某月平均气温小于或等于-10时,低温胁迫因子取0
        IF tmp_mon_cellvalue LE -10.0 THEN BEGIN
          e_t1[yindex,xindex] = 0
        ENDIF ELSE BEGIN
          ;低温胁迫因子
          e_t1[yindex,xindex] = 0.8 + (0.02 * tmp_opt_cellvalue) - (0.0005 * (tmp_opt_cellvalue^2))          
        ENDELSE       
      ENDELSE
           
    ENDFOR
  ENDFOR
 
END